Modeling,  Fabrication,  and  Electrical  Testing  of  Metal- 

Insulator-Metal  Diode 


by  Dr.  Terrance  O’Regan,  Matthew  Chin,  Cheng  Tan, 
and  Dr.  Anthony  Birdwell 


ARL-TN-0464 


December  2011 


Approved  for  public  release;  distribution  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position 
unless  so  designated  by  other  authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an  official  endorsement  or 
approval  of  the  use  thereof. 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


Army  Research  Laboratory 

Adelphi,MD  20783-1 197 


ARL-TN-0464 


December  2011 


Modeling,  Fabrication,  and  Electrical  Testing  of  Metal- 

Insulator-Metal  Diode 


Dr.  Terrance  O’Regan,  Matthew  Chin,  Cheng  Tan, 
and  Dr.  Anthony  Birdwell 
Sensors  and  Electron  Devices  Directorate,  ARL 


Approved  for  public  release;  distribution  unlimited. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the 
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Introduction 


Metal-insulator-metal  (MIM)  diodes  that  rely  on  quantum-mechanical  tunneling  for 
current  generation  are  one  of  the  few  structures  capable  of  rectifying  radiation  absorbed  by 
a  nano-antenna  designed  for  the  terahertz  frequency  range  (visible  and  infrared 
wavelengths)  (1,  2).  For  this  application,  diodes  with  current- volt  age  (I-V)  characteristics 
that  show  high  asymmetry,  strong  nonlinearity,  fast  responsivity,  low  hysteresis,  and  low 
turn-on  voltage  are  desired  (5).  We  show  that  these  properties  depend  primarily  upon  the 
materials  that  make  up  the  MIM  stack. 

We  model  the  tunneling  and  thermionic  currents,  and  compare  to  experimental  data 
measured  on  MIM  diodes  that  we  fabricated.  Our  model  gives  insight  into  the  quality  and 
variability  of  the  fabricated  MIM  structures.  The  model  confirms  that  the  MIM  diode  must 
be  dominated  by  quantum  mechanical  tunneling  and  be  strongly  asymmetric  to  rectify 
terahertz  radiation.  We  use  our  model  to  investigate  the  design  space  that  meets  these 
requirements. 

The  outline  of  this  report  is  as  follows.  First,  we  discuss  the  potential  energy,  an 
approximation  for  the  image  potential,  and  the  total  potential  energy.  We  show  how  the 
image  potential  lowers  the  total  potential  barrier  and  its  dependence  on  dielectric  constant 
of  the  insulator.  Then,  we  discuss  the  thermionic  emission-limited  current  as  a  function  of 
maximum  barrier  height  and  temperature.  Next,  we  explain  tunneling  limited  current  and 
its  dependence  on  average  barrier  height  and  tunneling  distance,  and  explore  thermionic 
emission  and  tunneling  current  as  a  function  of  insulator  thickness  and  metal-insulator 
barrier  heights.  Finally,  we  describe  the  MIM  fabrication  process  and  compare  our  model 
to  experimental  data. 


MIM  Model 


2.1  Potential  Energy  and  Image  Potential 

The  magnitude  of  the  intrinsic  electric  field  between  two  electrodes  separated  by  an 
insulating  film  is 

jp  _  (^ml  —  &m2) 

es 

where  <Fmi  and  <f>,ra2  are  the  metal  workfunctions,  s  is  the  insulator  thickness,  and  e  is  the 
unit  of  electron  charge.  The  barrier  heights  at  the  metal-insulator  interfaces  are  related  as 


<^2  =  ^i  +  eFiS, 
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where  <p>\  and  Lp2  are  the  metal-insulator  barrier  heights.  It  is  straightforward  to  show  that 
the  potential  energy  between  the  electrodes  can  be  then  written  as 

<p(x)  =  <p!  + A$(x/s),  (3) 

where  A<f>  =  <3>mi  —  <f>m2  and  x  is  the  position  coordinate  where  we  have  defined  the  diode 
to  be  oriented  along  the  x-axis.  If  we  assume  that  clectrode-1  is  grounded,  and  a  voltage  is 
applied  to  elcctrode-2,  the  potential  energy  between  the  electrodes  is  now 

(p(x)  —  ipi  +  (A<f>  -  eV)(x/s),  (4) 

where  V  is  the  applied  bias.  To  hnd  the  barrier  height  between  an  arbitrary  metal  and 
insulator,  we  need  the  workfunction  of  the  metal,  <Fm,  and  the  electron  affinity  of  the 
insulator,  Xins.  The  barrier  height  then  follows  as 

P  =  Xins-  (5) 

We  show  later  that  the  barrier  heights  are  critical  in  determining  the  current  through  an 
MIM  diode. 


The  image  potential  between  two  electrodes  separated  by  an  insulator  can  be 
approximated  as  (4,  5) 


Tim  (x )  — 


—0.288s 


Kx  (s  —  x)  ’ 

where  s  and  x  are  in  units  of  nm,  and  K  is  the  high-frequency  dielectric  constant  of  the 
insulator. 


(6) 


The  total  potential  energy  in  the  insulator  is  then 

htot  =  y{x)  +  Vim(x), 


or  more  explicitly 

htot  =  P>i  +  (A<&  -  eV)(x/s)  +  —— - 

Kx  (s  —  x 


(7) 

(8) 


Figure  1  shows  the  potential  energy,  image  potential,  and  total  potential  energy  for  a 
symmetric  diode  with  ipi  =  (£2  =  1  eV,  s  =  10m,  K  =  4,  and  V  =  1  V.  If  we  denote  the 
Fermi  levels  in  elcctrode-1  and  clectrode-2  as  Epi  and  Ep2,  respectively,  the  applied  bias  is 
Ep  1  —  Ep2-  The  image  potential  is  negative  and  reduces  the  total  maximum  potential 
barrier.  The  image  potential  also  smooths  out  the  abrupt  corners  of  the  potential  energy. 

Figure  2  shows  the  potential  energy  and  total  potential  energy  for  the  same  symmetric 
diode  as  in  figure  1  but  for  varying  values  of  K .  As  K  increases,  the  area  under  the 
potential  energy  barrier  increases  and  the  potential  energy  barrier  tends  towards  (p.  Thus, 
for  insulators  with  smaller  dielectric  constants,  the  peak  barrier  and  average  barrier  heights 
decrease,  which  we  show  later  increases  the  thermal  emission  and  tunnel  currents. 
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Figure  1.  Potential  energy  at  1-V  forward  bias  for  a  symmetric  MIM  diode  with  f\  =  1  eV, 
if 2  =  1  eV,  s  =  10  nm,  K  =  4,  m  =  mo,  and  T  =  300  K. 


Figure  2.  Total  potential  energy  with  varying  dielectric  constant,  K,  values  for  the  symmetric  MIM 
diode  in  figure  1. 


2.2  Thermionic  Emission-limited  Current 

The  thermionic  emission-limited  current  through  an  MIM  diode  is  (5) 

J  =  AthT2e~^/kBT  (1  -  e“eV/kBT)  ,  (9) 
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where  ip\  is  the  maximum  barrier  height  above  the  negatively  biased  electrode,  T  is  the 
temperature,  and 

.Ath  =  47t  mekx/h3,  (10) 

where  m  is  the  electron  mass  (set  to  m0  throughout  this  report),  ks  is  Boltzmann’s 
constant,  and  h  is  Plank’s  constant. 


From  equation  9,  we  see  that  the  thermionic  emission-limited  current  is  exponentially 
dependent  on  (always  referenced  to  the  Fermi  level  of  the  negatively  biased  electrode). 
From  figure  3,  we  see  that  is  constant  (set  by  —  Xms)  when  neglecting  the  image 
potential,  ffowever,  when  including  the  image  potential,  ip[  decreases  with  increasing  bias. 
In  fact,  at  1-V  forward  bias,  is  nearly  0.2  eV  below  the  no  image  potential  case. 

Figure  4  shows  the  calculated  thermionic  emission-limited  current  through  a  symmetric 
MIM  diode  described  in  figure  1  with  and  without  the  image  potential  included.  The 
conclusion  here  is  that  the  image  potential  must  be  included  if  one  hopes  to  accurately 
model  the  thermionic  emission  limited  current  through  an  MIM  diode. 
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Figure  3. 


Maximum  energy  barrier  above  the  negatively  biased  Fermi  level  for  the  forward  biased 
MIM  diode  in  figure  1. 
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Figure  4.  Thermionic  emission-limited  current  through  the  symmetric  MIM  diode  in  figure  1. 

2.3  Tunneling-limited  Current 

The  tunnel-limited  current  through  an  MIM  diode  is  (4) 


where 


J  =  J0 1  <p  exp  [~Atunip  2 
—  (</?  +  eV)  exp  -Aun  (<P  +  eV)~- 

f=> 


(n) 


(12) 


2nh(PAs)2> 

and  the  mean  barrier  height,  (p,  above  the  Fermi  level  of  the  negatively  biased  electrode  is 


1  C2 

<P  = 


(13) 


where  Si  and  s 2  are  the  positions  where  the  potential  energy  barrier  intersects  the  Fermi 
level  of  the  negatively  biased  electrode  and  As  =  S2  —  si  is  the  tunneling  distance.  The 
Atun  factor  is 

Aun  =  (47r/?As//i)  (2m)  2  ,  (14) 

where  f3  is  a  correction  factor  given  in  the  appendix  of  reference  4  (we  have  set  /3  —  1  in 
this  report). 


In  figure  4,  thermionic  emission-limited  current  dominates  for  s  =  10  nm.  Now  let  us 
explore  a  symmetric  MIM  diode  with  the  same  parameters  as  figure  1  but  with  s  =  5  nm. 
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Figure  5  shows  the  potential  energy  and  total  potential  energy  with  the  image  potential 
included.  Compared  to  figure  1  (s  =  15  nm),  the  total  potential  barrier  is  significantly 
reduced.  Figure  6  shows  that  the  maximum  potential  barrier  height  is  higher  than  the 
average  potential  barrier  height  for  each  applied  bias.  Comparing  figure  6  to  figure  2,  we 
see  that  the  maximum  barrier  height  is  reduced  as  s  is  reduced  due  to  the  larger 
contribution  of  the  image  potential.  Figure  7  shows  the  absolute  value  of  tunnel-limited, 
thermal-emission  limited,  and  total  current  versus  applied  bias.  At  low  bias,  the  thermal 
emission-limited  current  dominates.  The  tunneling-limited  current  increases  rapidly  with 
bias  and  crosses  over  the  thermal  current  at  higher  bias. 


Figure  5.  Potential  energy  at  1-V  forward  bias  for  a  symmetric  MIM  diode  with  ipi  =  1  eV, 
(^2  =  1  eV,  s  =  5  nm,  K  =  4,  m  =  too,  and  T  =  300  K. 
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Figure  6.  Maximum  and  average  barrier  height  vs.  applied  bias  for  the  symmetric  MIM  diode  in 
figure  5. 


Figure  7.  Absolute  value  of  tunnel-limited,  thermal  emission-limited,  and  total  currents  vs.  applied 
bias  for  the  symmetric  MIM  diode  in  figure  5. 

Next  we  consider  a  highly  asymmetric  diode  with  one  nearly  ohmic  metal-insulator 
junction.  In  this  case,  we  assume  =  0.2  eV,  </?2  =  1-2  eV,  s  =  5  nm,  K  =  4,  and 
T  =  300  K.  Figure  8  shows  the  absolute  value  of  the  current-density  components.  There 
are  three  important  characteristics  in  figure  8  that  are  important.  First,  the  tunneling 
current  is  not  nearly  as  asymmetric  as  the  thermal  current;  second,  the  tunnel  current 
dominates  the  thermal  current  for  all  negative  voltages  but  only  for  small  positive  voltages; 
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and  third,  the  tunneling-limited  current  acts  as  a  diode  in  the  reverse  polarity  compared  to 
the  thermal-limited  current.  This  second  point  may  not  pose  a  problem,  as  the  diode  will 
function  at  lower  AC  (terahertz)  bias.  However,  the  tunnel  current  must  be  made  more 
asymmetric  for  the  MIM  structure  to  act  as  a  diode.  The  third  point,  sometimes  referred 
to  as  “rectification  reversal”  is  due  to  the  switching  between  Fowler-Nordheim  and  direct 
tunneling.  Rectification  reversal  can  be  explained  by  figure  10,  where  the  average  barrier 
height  and  tunneling  distance  have  offset  each  other,  which  leads  to  a  reversed  polarity 
rectification,  and  in  general,  a  more  symmetric  I-V  characteristic.  Figure  9  show  the 
current  density  for  the  same  parameters  in  figure  8,  but  with  =  2.2  eV.  In  this  case,  the 
tunneling-limited  current  clearly  dominates  (thermal  current  is  not  shown),  but  the  I-V 
characteristic  is  still  quite  symmetric. 


Figure  8.  Absolute  value  of  tunnel-limited,  thermal  emission-limited,  and  total  currents  vs.  applied 
bias  for  an  asymmetric  MIM  diode  with  ipi  =  0.2  eV,  ip2  —  1-2  eV,  s  =  5  nm,  K  =  4, 
m  =  mo,  and  T  =  300  K. 
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Figure  9.  Absolute  value  of  tunnel-limited,  thermal  emission-limited,  and  total  currents  vs.  applied 
bias  for  an  asymmetric  MIM  diode  with  ipi  =0.2  eV,  </?2  =  2.2  eV,  s  =5  nm,  K=  4,  m  =  too, 
and  T  =  300  K. 
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Figure  10.  Maximum  barrier  height  and  tunneling  distance  vs.  applied  bias  for  the  asymmetric  MIM 
diode  in  figure  9. 


3.  MIM  Diode  Fabrication 


MIM  diodes  were  designed  and  fabricated  for  the  purpose  of  validating  the  model  discussed 
in  the  previous  sections.  Simple  column  structures  composed  of  a  blanket  bottom  metal 
layer,  blanket  insulator  layer,  and  a  square-patterned  top  metal  layer  were  chosen  due  to 
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their  ease  of  fabrication.  An  image  of  the  photolithographic  mask  design  for  the  column 
structures  is  shown  in  figure  11. 


Figure  11.  Simple  square  pillar  structures  were  used  as  MIM  diodes  based  on  niobium  (Nb) /niobium 
pentoxide  (Nb205)/x  material  stacks,  where  x  designates  copper  (Cu),  silver  (Ag),  and  gold 
(Au)  for  different  device  sets.  The  square  pillars  had  dimensions  of  60  fa n  x  60  /im,  and 
were  drawn  in  AutoCAD. 

The  Nb/Nb2C>5  material  system  was  chosen  due  to  this  combination’s  reported  low  barrier 
height  (6,  7),  which  provides  a  low  turn-on  voltage  in  a  MIM  diode.  Previously  reported 
measurements  have  specified  the  barrier  height  of  the  Nb/N^Os  junction  to  be  between 
133  to  300  meV  (6,  7).  The  top  metals  were  chosen  for  their  relatively  high  work  function 
difference  with  respect  to  Nb.  A  high  work  function  difference  between  the  two  metal 
layers  is  considered  the  primary  mechanism  for  high  asymmetry  and  large  nonlinearity  in 
MIM  diode  I-V  characteristics  (3,  8).  It  has  already  been  reported  that  Cu,  Ag,  and  Au 
have  relatively  good  asymmetric  I-V  characteristics  when  acting  as  the  second  metal  layer 
in  a  Nb/M^Os/x  MIM  diode,  where  x  represents  one  of  the  three  metals  listed  (7). 
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A  blanket  80-nm-thick  Nb  metal  film  was  deposited  via  a  room-temperature  sputter 
process  on  top  of  a  500- /un-t  hick  silicon  (Si)  substrate.  Next,  15-nm  of  Nb20s  was  grown 
via  an  electrochemical  (anodic)  oxidation  process  on  top  of  the  Nb  metal  film,  with 
approximately  5  nm  of  native  Nb20s  already  present  due  to  the  exposure  to  ambient  air 
after  the  Nb  film  deposition.  The  surface  roughness  of  the  as-deposited  Nb  film  was 
measured  to  be  0.9  nm  (average,  root  mean  square  [RMS]);  and  after  the  Nb20s  films  were 
electrochemically  grown,  the  surface  roughness  was  measured  to  be  1.1  nm  (average,  RMS). 
The  Nb205  thin  films  were  provided  by  collaborators  at  the  National  Renewable  Energy 
Laboratory  (NREL)  in  Golden,  CO.  A  WITec  Alpha  300  A  atomic  force  microscope 
(AFM)  was  used  to  characterize  the  surface  roughness  of  the  Nb  and  Nb20s  layers 
(measurements  performed  at  NREL).  The  thickness  of  the  Nb  metal  film  was  determined 
via  profilometry,  while  the  thickness  of  the  Nb20s  film  was  determined  using  cllipsometry, 
transmission  electron  microscopy  (TEM)  cross-sections,  and  x-ray  reflectivity. 

The  top  metals  were  photolithographically  patterned  using  a  lift-off  process.  For  model 
validation  purposes,  the  Au,  Cu,  and  Ag  top  metals  were  deposited  on  separate  sample 
sets.  The  geometric  structures  were  photo-patterned  using  the  AZ5214E  photoresist  and  a 
Karl  Suss  MA6  photoaligner.  Au  and  Ag  were  deposited  using  a  CHA  electron-beam 
evaporation  system,  with  a  deposition  rate  of  2  A/s.  Cu  was  deposited  using  a  Kurt  Lesker 
PVD75  thermal  evaporation  system,  with  a  deposition  rate  of  3  A/s.  For  each  of  the  three 
top  metals,  approximately  1,600  A  was  deposited  on  top  of  the  Nb20s  layer.  The  thickness 
of  the  top  metal  films  was  determined  through  the  use  of  profilometry.  Figure  12  provides 
an  example  of  an  optical  microscopic  view  (from  the  top)  of  a  device  set  fabricated  using 
the  process  described  in  this  section,  while  figure  13  provides  a  cartoon  image  of  a 
completed  MIM  diode  stack. 


4.  Comparison  to  Experimental  Measurements 


The  fabricated  devices  were  measured  with  a  Keithley  4200  Semiconductor 
Characterization  System  connected  to  a  Desert  Cryogenic  Probestation.  The  experimental 
measurements  were  done  at  a  temperature  of  298  K  and  vacuum  pressures  below  1  x  10-4 
torr.  The  electrical  characterization  involved  a  dual  voltage  sweep  from  —1  to  1  V  and 
back  in  voltage  steps  of  1  x  10-2  V  on  the  top  electrode  while  holding  the  bottom  electrode 
constant  relative  to  ground.  A  2.0-s  hold  time  was  used  prior  to  the  first  current 
measurement  on  each  sweep  to  allow  devices  to  settle. 
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Figure  12.  The  image  shown  is  an  example  of  MIM  diodes  fabricated  using  the  process  described  in 
this  section. 

The  resulting  current  was  then  extracted  and  used  for  fitting.  During  the  fitting  process, 
parameters  </?i  and  (f2  were  modified  in  order  to  produce  a  line  of  best  fit.  Parameters  T 
and  s  were  set  to  values  of  298  K  and  18  nm  (to  account  for  15  nm  of  grown  oxide  and 
5  nm  of  low-quality  native  oxide),  respectively.  On  the  other  hand,  K  was  experimentally 
determined  for  each  individual  I-V  measurement  from  the  slope  of  the  In  a  vs.  Ek b  *  T 
graph  assuming  electrode-limited  conduction,  where  a  is  the  conductivity  and  E  is  the 
electric  field  ( 9 ).  If  the  value  of  K  was  experimentally  determined  to  be  less  than  1,  its 
value  was  set  to  1  in  the  fitting  process. 
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Ml M TUNNEL  DIODE 


Figure  13.  In  this  cartoon,  a  top  and  cut-away  view  of  the  MIM  tunnel  diodes  fabricated  are  provided 
with  the  material  stack  thicknesses  and  dimensions  of  the  pillar  structure. 

The  results  of  the  three  different  devices  composed  of  Nb/Nb^Os/x,  where  x  can  be  Cu, 
Ag,  or  Au  are  shown  in  figures  14,  15,  and  16.  The  values  of  the  modeling  parameters  used 
are  given  in  table  1.  Nb  is  the  bottom  electrode  in  all  the  experiments  performed. 


Figure  14.  Modeled  and  experimental  current  characteristics  for  an  Nb/Nb205/Cu  device. 
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Figure  15.  Modeled  and  experimental  current  characteristics  for  an  Nb/Nl^Os/Ag  device. 


Figure  16.  Modeled  and  experimental  current  characteristics  for  an  Nb/Nl^Os/Au  device. 
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Table  1.  Best-fit  MIM  modeling  parameters. 


Mi/I/Mj 

K 

<Pi  (eV) 

y2  (eV) 

Nb/Nb205/Cu 

1.36 

0.55 

0.7 

Nb/Nb205/Ag 

1.00 

0.77 

0.94 

Nb/Nb205/Au 

1.60 

0.66 

0.81 

From  the  fitting,  we  can  see  that  there  are  slight  differences  between  the  experiment  and 
theoretical  data,  though  the  model  gives  a  good  match  to  the  experimental  behaviors  in 
general.  The  differences  are  primarily  detected  in  the  positive  bias  response,  which  is 
mostly  determined  by  the  Nb/N^Os  barrier  interface  tpi.  This  suggests  that  there  are 
parameters  or  non-idealities  in  the  device  unaccounted  for  by  the  modeling.  Ideas  for  such 
non-idealities  include  trapped  charges  in  the  oxide  and  noise  introduced  from  fabrication 
steps. 


5.  Conclusions  and  Future  Work 


In  the  present  work,  we  fabricated  MIM  diodes,  electrically  tested  them,  and  used 
modeling  to  extract  high-frequency  dielectric  constants  and  metal-insulator  barrier  heights. 
In  addition,  the  modeling  gave  insight  into  the  design  space  necessary  to  produce 
asymmetric  tunneling  I-V  characteristics.  Specifically,  thermionic  current  was  shown  to  be 
limited  by  large  barrier  heights  while  tunneling  current  was  shown  to  dominate  when  using 
thin  (~  5  nm)  dielectrics.  We  showed  that  the  metal  workfunctions  must  also  be  highly 
asymmetric  in  order  to  produce  asymmetric  I-V  curves.  Asymmetric  I-V  characteristics 
dominated  by  quantum-mechanical  tunneling  were  shown  to  be  difficult  to  achieve  because 
of  the  rectification  reversal  effect. 

In  future  work,  we  plan  on  verifying  the  electrically  extracted  barrier  heights  through 
alternative  microanalysis  techniques  and  by  investigating  other  material  systems.  Surface 
photovoltage  spectroscopy,  photoemission,  and  thermal  emission  measurement  methods  will 
be  used  to  confirm  the  logarithmic  electrical-conductivity  measurements  collected  in  this 
report.  The  results  from  these  methods  will  be  inserted  into  the  current  model  to  increase 
accuracy  and  provide  additional  insight  into  the  electrical  properties  of  MIM  diodes. 

Additionally,  future  MIM  diodes  will  be  fabricated  using  thinner  dielectric  films  with 
thicknesses  ranging  from  2  to  5  nm.  This  will  allow  us  to  investigate  the  contributions  of 
both  tunneling  and  thermionic  emission.  Low  temperature  and  graded  temperature 
measurements  will  provide  insight  to  the  contribution  ratio  between  the  two  emission  types 
in  MIM  diodes.  Temperature  measurements  and  calculations  will  provide  the  ability  to 
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further  extract  metal-insulator  barrier  heights.  Different  metal  contact  combinations  and 
other  metals  will  be  reviewed  to  determine  possible  material  systems  that  provide  even 
higher  asymmetric  MIM  diode  I-V  characteristics  than  current  MIM  diode  devices.  Besides 
the  Nb,  Cu,  Ag,  and  Au  reviewed  in  this  report,  platinum  (Pt),  tantalum  (Ta),  nickel  (Ni), 
and  chromium  (Cr)  may  also  be  investigated. 
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